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ContextThe increasing data set of precise observations of very energetic 
and collimated jets, with black hole (BH) as putative central engine, at 
different astrophysical scales and in various environments, should soon 
permit to discriminate and classify current theoritical models able to de- 
scribe the jets formation. 

Aims Constructing a purely gravitational theoretical model of perfectly 
collimated jets of high energy particles in the ideal case where the central 
engine is a Kerr BH of mass M and angular momentum by unit of mass 
a. 

Methods Studying in Weyl coordinates (p, z) the unbound Kerr 2D- 
geodesics which are asymptotes to straight lines parallel to the z axis of 
equations 



of which existence was recently demonstrated (Gariel et al.,2010). On 
these geodesies, flow test particles of energy E, with a Carter constant Q 
and (necessarily) an angular momentum L z — 0. 

Results We express the motion constants E and Q as functions of ri 
and r2, which are real roots of characteristics of the geodesies equations 
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Abstract 




system. In the special case of a double root r\ — r2 = Y, and, as an 
example, fixing a = M/2, the Y parametrization of the constants E and 
pi displays the following properties: 1) When E — > oo it implies \ Q\ — > oo, 
but pi remains finite and tends to only two possible values, pi/ pe = 10.241 
and 0.69 with p e = a/M. 2) E steeply decreases from infinity to small 
values, while pi concomitantly varies little inside two narrow ranges of 
pi/p e : [10.241; 10.65] and [0.69; 0.67]. Thus, the jet has a radial struc- 
ture. Hence, the energy flux can be calculated. Furthermore, based on 
observed data of the jets powers, we can obtain the mean particles density, 
the particles flow, the speed and the Lorentz factor of the jets, for any 
charged or neutral test particle. Then, we numerically apply these results 
to electrons. By studying the characteristics, we discuss the domains of 
initial conditions for geodesies starting inside the ergosphere. All these 
results come from the Kerr spacetime structure, and enhances the Penrose 
process as a plausible origin for the high energy jets. 

keywords: Astrophysical jets; Kerr black hole geodesies; High energy cosmic 
rays 

1 Introduction 

Among the numerous observed outflows ejected from various astrophysical struc- 
tures at all scales [2H1 02] the longest, the most energetical and collimated jets 
emerge either very fugaciously and intensely from cosmological sources - like 
the long duration gamma ray bursts [T] -, or with lifetimes at our time scale, 
steadily or repeatedly, nearer, from some microquasars, or, at much larger time 
and space scales, from the active galaxies nuclei [35]. All of them are relativistic, 
and often ultra relativistic [35]. Most of these jets are believed to be powered 
by a central engine, being a neutron star or BH, fed by larger structures, like 
elliptical galaxies, giant companion stars and X-ray binaries [34j . The more 
and more numerous and precise |26j observations of these various jets will soon 
help to strongly constraint the various conjectured theoretical models. In this 
perspective, we are here interested to suggest a purely gravitational theoretical 
model of the formation of a highly energetic and collimated jet powered by a 
rotating BH. 

As a first approximation, assuming the system to be axisymmetrical and in 
stationary rotation, we can represent a jet as a set of test particles following 
Kerr's unbound geodesies focusing at infinity along the z axis. 

In this ideal framework, if we succeed to describe a perfectly collimated jet 
with high energies, the model will allow us to build far more realistic descrip- 
tions by taking the ambient medium, i.e. the matter (including a magnetic 
field), into account. However, these more realistic descriptions, although im- 
portant, will only produce marginal improvements concerning the origin of the 
jet formation. The essential focus of the phenomenon will remain the coherent 
set of parallel unbound geodesies of a Kerr BH combined with the source of 
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powering essentially being the gravitational field in its strongest manifestations, 
namely the BH. 

The generalized cylindrical, or Weyl, coordinates (p, z, 0), related to Boyer- 
Lindquist generalized spherical coordinates (r, 9, <f) by 

p= [(r- l) 2 - A} 1/2 sm9, z = (r-l)cosQ, (1) 

where 



are the most suitable for describing observable phenomena generated by ax- 
isymmetric structures. The existence of special unbound geodesies was recently 
demonstrated in this framework [14] . These geodesies stem from the ergosphere 
and, when z — > oo, they are asymptotically parallel to the z axis, with 

/ Q \ 1/2 

P = *= [6 + (3) 

for the asymptotes, depending on p e = a/M and on the two constants of motion, 
the Carter constant Q and the energy E, the third constant of motion, the z 
component of the angular momentum, L z being necessarily null. In the present 
paper, we show that only some of these geodesies, belonging to narrow ranges, 
can be followed by particles with high energies. 

The function R(r) [|Tj introduced in the expression of the Kerr timelike 
geodesies (test particle mass v^i = 1) plays a fundamental role in the anal- 
ysis of the jet collimation in the case of a model where the engine at the centre 
of the accretion disk is supposed to be a stationary rotating BH. This function 
is such that 

R 2 (r) = a 4 r 4 + a 3 r 3 + a 2 r 2 + air + a Q , (4) 
with (see (2) and (4-8) in [14]) 

ao = -a 2 Q, ai = 2{a 2 E 2 + Q), 

a 2 = a 2 (E 2 -l)-Q, a 3 = 2, a 4 = E 2 - 1, (5) 

where we put M = 1 and L z = 0, considering the special 2D-geodesics of ([3])- 
type. Hence, the BH spin a being fixed, — 1 < a < 1, we have two independent 
parameters left, Q and E, or, equivalently from ([3]), the position pi of the 
asymptote parallel to the z-axis and the energy E. 

The paper is organized as follows. In section 2, we obtain the expressions 
of the two motion constants, E and Q, as functions of two real roots of the 
characteristics equation R 2 (r) — 0. In section 3, we consider the special case 
of a double root Y , and we show that there exist only two narrow ranges of Y 
for which E can have high energy values. In section 4, we show that the two 
possible other roots are functions only of the two first ones. Then we look at 
some consequences on the admissible values of E and the corresponding ranges 
for the asymptotes p\. In section 5, we calculate the energy flux of the jet, and 
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based on the observational evaluation of the power of the jet we deduce the 
corresponding particles density, the particles flow, the mean velocity and the 
mean Lorentz factor of the jet. As an example, we give a numerical estimation 
of these quantities for electrons. In section 6, by studying the characteristics, we 
show that among the two previous possibilities found in section 3, there remains 
one admissible only. In section 7, as a conclusion, we discuss qualitatively some 
potential consequences of relaxing some restrictive assumptions made here on 
the possibilities offered by the Penrose process to obtain high energies with 
efficient jets formation. 



2 Conserved quantities as functions of two roots 

Let us consider the possible roots of the equation R 2 (r) = of the characteristics 
f = of the autonomous system of geodesies equations [S] , i.e. 

a 4 r 4 + a 3 r 3 + a 2 r 2 + air + a = 0. (6) 

The polynomial equation (j6)) has four roots, labeled 7"j with i = 1, 2, 3, 4, which 
can be a priori > or < or complex (contrarily to the r physical variable 



which is real defined in the range 



). The two equations R 2 {r\) = 



1 + VA, oo 

and R 2 (r 2 ) = are linear in Q and in E 2 — 1. Solving the linear system of these 
two equations yields the two parameters as functions of the roots n and r 2 , 

Q = {a 4 + a 2 [n(n - 2) + r 2 (r 2 - 2)] (7) 

+rlr 2 }, 

E 2 -l = -l{a* + a 2 (r 2 +r 2 2 ) (8) 



D 

-nr 2 [ri(r a -2)-2r a ]} > 



with 



D = a i (2 + r 1 + r 2 ) + a 2 [rf + r\r 2 + r x r 2 {r 2 -4) + rf\ 

+rir 2 [(rf + nr 2 )(r 2 - 2) - 2r|] . (9) 

In the third possible equation, R 2 (r^) = 0, the parameters Q and E 2 — 1 
can be replaced by ([7]) and ©, leading to a relation between r^ and r\ and r 2 
allowing, in principle, to determine the values of r 3 as functions of t\ and r 2 
only, with a being fixed. The fourth possible equation, i? 2 (r 4 ) = 0, will not 
bring any new result because the roots r 3 and r 4 are the same. 

In (0) and (JSj , it is worth noting the symmetric role of r± and r 2 , and that Q 
and E 2 — 1 have the same denominator D, so that if, and only if, it cancels, we 
have E — > oo and \ Q\ oo, whereas pi, depending only on their ratio (see ©), 
tends towards a finite value. From ©,0 and (J5J we obtain the asymptotes 

Pi\ 2 _ (a 2 + r 2 )(a 2 - ri r 2 ){a 2 + r 2 ) 



p e J a 2 {a i + a 2 {r 2 1 +rl) + r 1 r 2 [r x (r 2 - 2) - 2r 2 ]}' 
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3 Roots r\ = T2 real 

For sake of simplification, we assume in this paper that there is a double real 
root fi = r2 = Y. Hence (|4|) can be rewritten as 



R 2 {r) = a 4 (r - Y) 2 (r 2 + Br + C), 
and ((7|) and (jHJl simplify to 

[a i + 2a 2 {Y -2)Y + Y 4 ]Y 2 



Q = 



i 4 (l + Y) + 2a 2 (Y - l)Y 2 + (Y - 3)y 4 ' 



E 1 



■ 2a 2 Y 2 + (Y — 4)F 3 



a 4 (l + Y) + 2a 2 (Y - l)Y 2 + (Y - 3)F 4 ' 



(11) 

(12) 
(13) 



When E — > oo, of course, \Q\ — > oo again. But, as already noted, their ratio 
tends towards a finite value, so that pi remains finite, 



2 \2 



(a 2 -Y 2 )(a 2 +Y 2 ) 
a 2 [a 4 + 2a 2 Y 2 + (Y - 4)F 3 ] 



(14) 



Choosing a "moderate" rotation of the BH, we fix a = M/2 (as considered 
by us in [H]). Hence, we can plot the functions E 2 — 1 = F(Y) and (pi/p e ) 2 = 
G(Y), as shown in figures Q] and [2j 
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Figure 1: Plot of E 2 — 1 = F(Y), where E is the energy of the test-particle, 
in function of the double root Y, evaluated from the relation (I13p . for a BH of 
mass M = 1 and of angular momentum by unit of mass a/M — 0.5. We can see 
the ranges of Y for which E 2 — 1 is positive, as expected for unbound geodesies, 
and the three values of Y for which E 2 — 1 tends to the positive infinity. 

Since F and G have to be simultaneously > 0, the only possible solutions 
correspond to the two ranges 

Y E [-0.5, Y 0a ], (15) 

Y € [Y Qb , 3.86971] , (16) 
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Figure 2: Plot of (pi/p e ) 2 = G(Y) as function of the double root Y, evaluated 
from the relation (|14p , for a BH of mass M = 1 and of angular momentum by 
unit of mass a/M = 0.5. As in figure Q] we can see the ranges of Y for which 
the function (pi/p e ) is positive. 

with Yoa and Y b the asymptotes of F(Y), for which i5 — > oo. We can nu- 
merically evaluate these asymptotes (namely the roots of the equation D = 0), 
yielding approximately Y" a — —0.241806 and Ygfc — 2.8832 (for more precise 
values see (f2"6l and ([23)1 respectively). 

Hence, there are two only possible values of p\ for which E — > oo, i.e. one 
for each range (fT5|) and (|T6| . For the limits of the two ranges ((T5|) and (fT6|) . 
y = Foa — £ and Y = Fob + £, when e — ► 0, we obtain the finite values, 
respectively, 

— ~ 0.693199, and — ~ 10.2411. (17) 

Pe Pe 

At the other extremity of the range (TIB"]) , i.e. for Y = 3.86971 where pi oo, 
we have E 2 — 1 = 0. And for the other extremity of the range (fT5")) . i.e. for 
Y = —0.5 where E 2 — 1 = 2, we have p\ = 0. 

The figure [3] summarizes these results. 

4 Roots r3 and r4 

Identifying (TTTj) with (Q} rewritten with the parameters E 12 — 1 = i^(Y) and 
(pi/ p e ) 2 — G(Y), without the explicit form of these functions of Y, given by 
(fi"3]l and (fT4| . yields the four relations, 

s-2y = -|, g + ^ = 2(sr 2 - 2yc), 

_r F 

1-G=16CY 2 , 2-G = A(C-2YB + Y 2 ), (18) 
linear in 1/F, G, B and C. 
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Figure 3: Some precise key values of the functions F(Y) = E 2 — 1 and G(Y) = 
(pi/Pe) 2 showing the ranges of Y for which these functions are simultaneously 
positive (bold lines). The negative values of these functions are in dotted lines. 



After eliminating 1/F and G in (fl8|) . we obtain B(Y) and C(Y), namely 

2(4F 2 - 1)[Y(4Y - 1) + Y - 1] 



C 



( 4 y2 _ ^2 _ \QY 2 {4Y - 1) 

(4r 2 -l) 2 + 16F(r-l) 



(19) 
(20) 



4[(4F 2 - l) 2 - 16y 2 (4F - 1)] ' 
Hence, (fTTj) can be rewritten as 

R 2 = a 4 (r - Y) 2 {r 2 - Sr + P) = a 4 (r - Y~) 2 (r - r 3 )(r - r 4 ), (21) 
where r3 and r 4 are the remaining roots, in general distinct, and 

S = r 3 +U = -B, P = rzU = C, (22) 



or 



^3 = 

r 4 = 



B+{B 2 -ACf' 2 
B-(B 2 -4C) 1/2 



(23) 
(24) 



where B(Y) and C(Y) are given by (U]) and ([2D]). 

The curves ^(Y) and ^(F), are plotted in figures [4] and [5j and are real for 
some ranges of Y only. In particular, in the range (| 15[) for 1", r3 and r 4 are 
not real. To have the expression r 2 + Br + C in Eq. (fTTj) real, where i? and 
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C are real, r% and r 4 have to be complex conjugated, i.e. r% = z = B\ + iC± 
and r 4 = z. Hence, the sign of the expression r 2 + Br + C = (r + Bi) 2 + C\ is 
always positive, and P = C = Bf + C\ > and S = -B = -2B X < 0. 

r 3 
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-200 

Figure 4: Plot of the root r% in function of the double root Y, evaluated from 
the relation (l23l) . for a BH of mass M — 1 and of angular momentum by unit 
of mass a/M = 0.5. 
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Figure 5: Plot of the root r 4 in function of the double root Y, evaluated from 
the relation ([24)) . for a BH of mass M = 1 and of angular momentum by unit 
of mass a/M = 0.5. 

In the range (|16|) for V, the two roots and r 4 are real, P — C is negative 
(which means two roots of opposite signs) and B = —S is positive. The portions 
of curves r 3 (Y) and r 4 (F) can be plotted on the range (fi"6|) . The most precise 
value we can numerically obtain for the left limit (where, in principle, E — > oo) 
of the range (fT6|) , is 

F 06 ~ 2.8832177419263523927462568785232847, (25) 
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allowing us to reach the maximal value E ~ 1.1 * 10 32 "only". Then, the corre- 
sponding real values r 3 = -5.898835521038341 and r 4 = 0.13240003718564175 
are obtained. 

Also, it is worth observing that E is steeply decreasing, for a weak variation 
e of Y from Yo6 (e > 0) or from Yo a (e < 0), while p\ is weakly increasing 
for this same small interval of Y , For example, when Y goes from Y~ob to 
2.922, the energy E is steeply decreasing from 10 30 to 3, while the position 
of the asymptote p\/ p e increases by a small amount from 10.24 to 10.68, which 
means a big concentration of the most energetic part (the " spine" ) of the beam 
immediately near, at the right hand side, of pi/p e — 10.24. At its left hand 
side, there is no more beam produced. 

Likewise, for the range (fl~5|) . the energy E is very steeply decreasing from oo 
to 6. The highest value of E (E max ~ 2 * 10 30 ) is obtained for the most precise 
value numerically obtainable for the right limit Y"o a , namely 



The corresponding asymptote pi/ p e , inside the ergosphere, decreases very slightly 
from 0.6932 to 0.6764. Here the jet is yet more concentrated just at the left of 
0.6932. While beyond its right side, there is no possible beam. 

As a result, our model predicts a radial structure of the jet, with a precise 
profile for its energy (or speed) distribution of the particles. A radial morphology 
has already been suggested from observations [15l [21] . 

5 Flux, particles density and Lorentz factor 

In the region z 3> pi of the jet, where the beam is quasi-parallel to the z axis, 
the vector density of (total) energy current is purely convective, je = Pev, 
where pe — nE, being E the mean energy of a particle and n is the number 
density, v — ve z its velocity. The power of the jet or, equivalently, the energy 
flux across a crown, cross-section of surface S, included between two radii p±Q 
and pii, is 



Y 0a ~ -0.241805810271953623405344008301523644957. 



(26) 




(27) 



where dS = dSe z and we supposed an homogeneous jet, n = constant. The 
speed of each particle is (see equation (49) of |14j). 




its kinetic energy is defined by 



E c = (E 2 -c i S 1 ) 1 ^, (29) 
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and the total energy E is linked to the reduced function E(Y), given in (TTBl . by 

E = -rw- ( 3 °) 

Hence we have 

vE = c 3 (E 2 - l) 1/2 ^. (31) 

Likewise, the " true" length p\ in function of the reduced coordinate p\ , itself 
function of Y by the equation (fl"4"|) . is 

* = (32) 



M 



Whence, 

^ = -r 



= 7rnM 2 p2 c 3^ [ Y \e 2 - !)i/2 d (Pi/Pe) 2 (33) 



in 



where the adimensional functions -E'(y) and (/3i/p e ) 2 (y) are given by the equa- 
tions (fT3"|) and The relation (|3"3")l can be used to calculate, for exam- 
ple, the power of the beam of the particles ejected along the geodesies with 
asymptotes included between p w / ' p e = 0.69319914385 and pn/p e — 0.676252, 
corresponding, as seen in section 4, to energies varying from the infinity (in 
fact, E - 2 * 10 32 ) to 5.8, and with Y range (HJ from Y Qa ~ -0.241806 to 
Y\ = —0.250. However this calculation would suppose knowing the particles 
density n inside the jet. Conversely, if the power (which is an observational 
data more easily attainable by other independent methods [3S1 [T31 [23 122] ) , is 
known, then from our model we can deduce the density n. 

Further assumptions have to be made about the ejected particles and the BH 
parameters. Let us consider a supermassive BH with mass M = 10 9 M Q (with 
the solar mass Mq = 3km) and angular momentum by unit mass a = M/2 , and 
electrons with restmass y/6i — 0.51 lMeV. Let us suppose that the jet power is 
V = 10 47 erg/s [36l [31] ■ We find with this data (the calculation can be made 
numerically from a series expansion of the integrant near Yo a , which is justified 
because \Y -Y 0a \< |Yi - Y 0a \ ~ 8.2 * lO" 3 < 1), 



Y \ E j2_ 1) i/2d{ P i/Pe)\ Y 



'Yoa 

Whence the density is deduced 

10 47 



dY 



0.264939. (34) 



3.29 * 10 28 



2.18 * 10 20 electrons/m . (35) 
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The flow of particles across the same crown of surface S is 



dN_ 
~dt 



mr J vdp 1 

27 on 

n *3.14* — * 10 32 * J, 



and we obtain, 



J = 



2.30373* 10" 2 , 



1/2 



dY 



dY 



and finally, 



^ = 1.06 3 52 * 10 52 clcctrons/s. 



From these results the mean energy by particle can be deduced, 

-l 



< E >= 
= 0.511* 



~dl \~dt 
0.264939 



2.30373* 10- 2 
which yields the (mean) Lorentz factor 



5.87673 MeV/electron, 



0.264939 



c 2 V5i J 2.30373* 10- 2 
11.5005, 



and the mean velocity 



c 1 



< E >< 



1/2 



0.996212 *c, 



(36) 



(37) 
(38) 



(39) 



(40) 



(41) 



corresponding to an average ultra relativistic jet. 

Now, we also can consider a narrow, more energetic, part of the jet only. As 
examples, let us consider the part included between p w and pi2 = 0.69319 * p e 
(or P13 = 0.693199 * p e ), corresponding to Y"oa and Y 2 = —0.24181 (or Y 3 = 
-0.241805810272), respectively, i.e. to thickness Sp 2 = pn - pio = -8.68546* 
10~ 6 * p e , for Y 2 - Y 0a = -4.18973 * 10~ 6 (or 5p 3 = p 13 - p w = -9.61453 * 
10~ 14 * p e , for Y 3 — Y) a = —4.63518 * 10 -14 ), and hence to ratio of flux surfaces 
S2/S = (fi\ 2 - P 2 )/(P 2 ! - Plo) = 5.18849 * 10- 4 (or S 3 /S = (p\ 3 - p 2 )/(p 2 1 - 
p 2 w ) = 5.74535* 10- 12 ). 
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Furthermore, let us (reasonably) assume that the power by unit surface 
crossed by the jet (i.e. the energetic flux) is constant: 



V_ = Vi L 
S Si 




(42) 



Hence, the part of the jet with the thickness Sp 2 (or Sp 3 ), will have a kinetic 
power V 2 = {S 2 /S)V = 5.18849* 10 43 erg/s (or V 3 = {S 3 /S)V = 5.74534* 10 35 
erg/s), a particle density n 2 = {S 2 / S){I / I 2 )n = 4.97047* 10 18 electrons/m 3 (or 
n 3 = (S 3 /S)(I/h)n = 5.23365*10 14 electrons/m 3 ), a flow of particles dN 2 /dt = 
1.26855*10 47 electrons/s (ovdN 3 /dt = 1.48 031*10 35 electrons/s), a mean energy 
< E 2 >= Vdi(h/J 2 )c 2 = 255.631 MeV (or < E 3 >= y/h{h/J 3 )c 2 = 2.42573 
TeV), a mean velocity v rn2 = 0.99999800204c (or v m3 = 0.9999999999999778c) 
and a Lorentz factor T 2 = I 2 /J 2 = 500.257 (or F 3 = I 3 /J 3 = 4.74703 * 10 6 ). 
Let us precise that I 2 = -6.02378 * 10~ 3 and J 2 = -1.20414 * 10~ 5 (or 7 3 = 
-6.33485 * 10" 7 and J 3 = -1.33449 * 10~ 13 ) are the integrals (03) and ([37]) 
respectively, in which the upper limit becomes Y 2 (or Y 3 ) instead of Y\. So, 
we can see that, in its core, the jet is as more strongly UR as nearer of pro- 
Let us precise that Y 3 , and all the corresponding values with the index "3", 
corresponds to the value of Y the nearest of Y^ a for which we are yet able to 
numerically evaluate the integrals ^3 and J 3 (and so, the different corresponding 
quantities ). Let us note that the conditions inside the jet for the most narrow 
channel we are able to evaluate are comparable to conditions inside a terrestrial 
particles accelerator (e.g. LHC), namely < E 3 >~ 2 TeV inside a channel of 
thickness ~ 8p 3 ~ 14.4218 cm (but here for electrons, while in the LHC there 
are protons. We can evaluate jets of protons as well, which would give us yet 
higher energies: factor ~ 2 * 10 3 ). The results are summarized in the table 1. 

As second example, we give the same evaluations for the second part of 
the very energetic jet, namely near the value Y ob ~ 2.88 l|2"5"|) . The results are 
summarized in the table 2. We shall not comment more precisely these last 
results, because they cannot correspond to effective jets, as we shall show in the 
next section. 

6 Geodesies for high energy jets 

Now let us look at the possible geodesies framing a jet for values of the energy 
and of the corresponding asymptotes p\ given in section 5. The existence of 
admissible initial conditions for such geodesies are provided by the existence 
and positions of the 2D-characteristics of the system of geodesies equations (see 
(2) and (3) in 14 ). Each characteristics curve delimits two separated parts 
(regions) in the plane (r, 9) as predicted by the theory [30] often applied to the 
so-called "qualitative analysis". Each part contains a set of geodesies, which 
can never cross the characteristics towards the other part. In Boyer-Lindquist 
coordinates, the characteristics are defined by the equations 



r = 0, 9 = 0, 



(43) 
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Yi 


Y 2 


Y 3 


Y 


-0.250 


-0.24181 


-0.241805810272 


Y - Y 0a 


-8.19419* 10" 3 


-4.18973* 10"° 


-4.63518* 10" 14 


E 


5.81378 


250.13 


2.37778* lO 5 


Spi/Pe 


-0.0169469 


-8.68546 * 10 _e 


-9.61453* 10" 14 


Spi (m) 


2.54204* 10 1U 


1.30282* 10 7 


14.4218* 10~ 2 


P (erg/s) 


10 47 


5.18849* 10 43 


5.74534* 10 a5 


I 


-0.264939 


-6.02378* 10~ a 


-6.33485* 10- y 


n (e~/m 3 ) 


-2.18711 * 10 2U 


-4.97047* 10 1S 


5.23365* 10 14 


J 


-2.30373* 10~ 2 


-1.20414* 10" 5 


-1.33449* 10" 13 


dN/dt (e~/s) 


1.06352* 10 52 


1.26855* 10 47 


1.48031 * 10 a5 


E m (Mev) 


5.87673 


255.631 


2.42573 * 10 6 


1 - v rn /c 


3.788* 10~ a 


1.99796* 10 B 


2.22* 10" 14 


r 


11.5005 


500.257 


4.74703 * 10 B 



Table 1: Values of the departures of the position from the asymptotes at high 
energies, the thickness of the jet, its density, its mean energy, its mean velocity 
and its (mean) Lorentz factor, as function of the parameter Y, for Y in the 
vicinity of Y 0a ([26]), where E -> oo and pi/p e -> 10.24106 (E max = 1.1 * 10 32 ). 
The integrals / and J are defined by (|34]l and ([37)) with the corresponding upper 
limit Yi (i = 1, 2, 3), the lower limit being Y 0a . 



which are equivalent, from equations (19) and (20) in [T3], to the system of 
algebraic equations, 

P = 0, 5 = 0, (44) 

the solutions of which, when they exist, are some values of fj of r (previously 
introduced in Sections (2) to (4)) and some values 9i of 9, respectively (which 
define circles and straight lines from the origin respectively) . 

In Weyl coordinates, p and z, these characteristics equations (|44"]) are equiv- 
alent to the equations (see (17) and (18) in [14]) 

. S(a 2 -A)z Sa 

P= apA ' Z = ^ (45) 



and 



respectively, where 



VM 

Each set of equations (|4"5)) and (|4*6"|) lead to, 

dz a 2 p 

~dp = ~ (c?-A)z 



Pa 3 p . Paz 
P= (a 2 -A)A > Z = ^' (46) 



A = (a + l) 2 a 2 + [i-) 2 z 2 . (47) 



(48) 
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Yi 


Y 2 


Y 3 


Y 4 


Y 


2.922 


2.833 


2.88321774193 


2.883217741926354 


Y - Y ab 


3.87823* 10" 2 


8.2258* 10" 5 


3.6473* 10" 12 


1.33227* 10" 15 


E 


3.027 


63.7027 


302 514 


1.50555* 10' 




0.411539 


8.49897* 10" 4 


3.76801* 10" 11 


1.77636* 10" 14 


5pi (m) 


6.17309* 10 11 


1.27485* 10 a 


56.5201 


2.66454* 10" 2 


P (erg/s) 


10 4V 


2.02457* 10 44 


8.97533 * 10' ib 


3.96648 * 10'" 


I 


49.1803 


2.21766 


4.668989* 10" 4 


7.2287* 10" e 


n (e~/m ,:l ) 


1.17177* 10 lt! 


5.26106* 10 ltj 


1.10759* 10 13 


3.16214* 10 11 


J 


8.33589 


1.74073.* 10" 2 


7.62157* 10" 1U 


3.63748* 10" 12 


dN/dt (e-/s) 


2.07028 * 10 b2 


1.94106* 10 4S 


1.78919* 10 3 ' 


2.43823 * 10'" 


E m (Mev) 


3.01481 


65.1005 


313100 


1.01536* 10 B 


1 — Vm/c 


1.4469* 10" 2 


3.1 * 10" b 


1.3318* 10" 12 


1.266* 10" 13 


r 


5.89983 


127.398 


612 719 


1.98701* 10 s 



Table 2: Same values as table 1 but in the vicinity of Y"o& (|25|) where E —> oo 
(E max = 2 * 10 30 ) and Pl /p e ~ 0.69. 



and 



dz 
dp 



{a 2 - A)z 



a 2 p 



(49) 



respectively, defining the two families of characteristics for the geodesies of type 
(21) in |14) in which we are interested, namely, ellipses (corresponding to r = 0) 
and hyperboles (corresponding to 8 = 0). Let us note that the product of the 
two derivatives (|48|) and (|49|) of these characteristics is — 1 (which confirms that 
they are orthogonal). 

The first ones, ellipses, exist when there are solutions r = fj = constant of 
(|48|) for V0, with r, > 1 + V~A or equivalently a — on = constant ( because 



1) with cti > VA. Then ggj) can be integrated yielding 



Ki, 



A 



(50) 



where K\ is an integration constant. Comparing (1501) with the equation (12) of 
[T4"] imposes K\ = 1. 

The second ones, hyperboles, exist when there exist solutions /i = /i; = 
constant of (|49t for any r, with /uf < 1. These are solutions of the equation 
S = 0, when L z = 0, with 



1 - 



a 2 (E 2 



(51) 



There are two possible cases, namely /x? = 1, then S = for any Q, or pf — 
— Q/[a 2 {E 2 — 1)] = 1 - (pi/pe) 2 < 1 , being positive defined only if Q <0, or 
equivalently, if p\ < p e . Then, for /xf = 1, we have z = r — 1 = a and p = 
for any r and (|49[) reduces to — > oo , and the characteristics being along the 
semi-axis z > VA. While for p 2 = -Q/[a 2 (E 2 - 1)] we have for ([49]) 



dz 
dp 



zp 



(52) 
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which can be integrated leading to 



K 2 



"I 1/2 



— I -A 

Hi 



(53) 



where K 2 is an integration constant. Comparing (1531) with (12) in |14) we have 
The expression (|53[) represents a family of hyperboles parametrized by 



Pi 



,V2 



yielding 



1 
A 



1-l^i 

Pe 



(54) 



(55) 



If the initial condition (IC) of a geodesies lies inside an ellipse of the type 
(|50)l . this geodesies cannot be an unbounded geodesies, and hence cannot go to 
infinity. So, the admissible IC have to satisfy the triple condition: i) being inside 
the ergosphere, in order to be possibly issued from a Penrose process; ii) being 
outside the larger elliptic characteristics, this means corresponding to the larger 
value of the roots r^; and iii) being above the higher hyperbolic characteristics, 
which corresponds to the higher values of the roots < 1. That restricts the 
admissible domain of IC. 

H 



An ellipse (f50|) . when it exists (i.e. when r, e 1 + \/ A,oo ), can intersect 

the ergosphere only if its semi-minor axis hi = [a 2 — A) 1 /' 2 is smaller than 
p e = a/M, i.e. if r, < 2. 

As example, let us take the special case of a double root Y studied in the 
precedent sections. 

a) The first admissible range that we found is Y E [io&>3.86] (see figure [3]). 

These roots, belonging to the domain of physical definition, re 1 + \[~A = 1.86, oo 

all correspond to the existence of elliptic characteristics. The smallest ellipse 
has as semi-minor axis b t = [{Y -l) 2 - A] 1 / 2 = [(1.88) 2 -0.75] x / 2 = (2.78844) 1 / 2 
along /5, and as semi-major axis a,i = \focf — Y — 1 = 1.88 along z, obtained 
for the smallest value Fob — 2.88, corresponding to pi/p e — 10.24. This ellipse 
contains the ergosphere, the limits of which being z max = \[A — 0.866025 and 
Pmax = 1/2. Hence, it is always impossible to have IC simultaneously inside the 
ergosphere and outside any ellipse. There is no possibility of unbound geodesies 
starting from the ergosphere in this first case. 

b) For the second admissible range we found that Y 6 [—0.5, Y 0a ] (see figure 
[3]). These roots do not belong to the domain of definition of the physical variable 
r, which means that there are never any corresponding elliptic characteristics. 
The only remaining possible limitation depends on the position of the hyperbolic 
characteristics (|55|) . The hyperbola intersects the z -axis at the point with 
coordinates p = and Zo = {A[l—(pip e ) 2 }} 1 / 2 and tends asymptotically towards 
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the straight line of equation p ~ ztan#i, with sin#i = p\/p e - The domain 
of possible IC is located between the z-axis, the limit of the ergosphere and 
above the hyperbola. For example, for Y = -0.241806, p 1 /p e = 0.693199, 
6»i = 21°, and z = 0.624185(< V%~ = 0.8660254037844386). We plot in figure[6] 
the geodesies which tends asymptotically towards the corresponding p = p\ = 
0.3466, for which the test particle has a very high (theoretically infinite) energy 
(for the calculations, we choose the value E = 10 6 and pi/p e = 0.693199). This 
plot corresponds to the IC p t = 2.8 * 10~ 6 and z, = 0.852086186870110 which 
are (just) inside the ergosphere at its top near the z-axis, i.e. near the event 
horizon. For the other limit, Y = —0.5, of the Y range, p\ = , E = v3> #i = 
and z Q = y/A = 0.866. 

P 

0.3466 
0.3465 
0.3464 




0.3463 



\ 20000 40000 60000 80000100000 

Figure 6: Plot of the geodesies p(z), for the parameters a = 0.5M and M = 1 
of the BH, the motion constants L z = 0, E = 10 6 , Q = [{pi/p e ) - 1] (E 2 - 
1) (a/M) 2 (with pi/ p e — 0.693199) of the test-particle, and for the initial con- 
ditions p t = 2.8 * 10~ 6 , Zi = 0.852086186870110, inside the ergosphere. We 
can observe that this geodesies is asymptotically parallel to the z-axis at the 
asymptote of equation p = p\ = 0.3466, at least until the altitude z ~ 8 * 10 4 . 
For a greater altitude, a greater precision on the initial conditions would be 
necessary, theoretically without any limitation. 



7 Discussion 

Taking as parameters the roots of a characteristics equation for unbound 2D- 
geodesics with L z = 0, we showed that the two remaining motion constants, E 
and Q, of a test particle following geodesies which asymptotically tends towards 
a parallel line to the z-axis, can be deduced as a function of these parameters. 
In the special case of a double root, and choosing the BH angular momentum 
by unit of mass a = M/2, restricted domains of p\ asymptotes corresponding to 
high energies are found. That means that the Kerr metric can generate powerful 
collimated jets in some precise regions only (and, as a consequence, from some 
precise regions of the ergosphere). Indeed, we obtained, in this special case, two 
only possible ranges of pi, namely pi £ [0.3382,0.3466] and p\ £ [5.12,5.34] 
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for E 6 [V6, oo[ and E g ]oo, v3] respectively (see the figure [3]). Then, for 
electrons, we calculated the particle density, the particle flow, the thickness of 
the jet and the Lorentz factor, exhibiting UR jets. The numerical evaluations 
were summarized in the tables 1 and 2. 

These results can be pertinent both for jets and for ultra high energy cosmic 
rays (UHECR). Indeed, such energies can be obtained with quasi unlimitedly 
values inside the ergosphere thanks to the Penrose process. That comes from 
the fact that this process implements three levels (orders of magnitude) of rest 
masses (or of energies): the BH mass (which could be designated as a "super- 
macroscopic", or astrophysical mass), the incident mass of the body going into 
the ergosphere along the equatorial plane from the accretion disk (which could 
be called a "macroscopic" mass); and the particle going out of the ergosphere 
asymptotically parallel to the z axis (which could be called a "microscopic" 
mass) . 

The "microscopic" outgoing particle, via the original Penrose process [27] , 
simply by its separation from the incident macroscopic body inside the ergo- 
sphere, can acquire a kinetic energy of order of magnitude of this ingoing mass, 
i.e. a "macroscopic" kinetic energy. 

The simplest example of such a process could be the ionization of an ingoing 
atom of hydrogen inside the ergosphere, with its nucleus irreversibly falling into 
the event horizon while the electron is ejected along a geodesic asymptotically 
to p\ with a kinetic energy increased of about 1836 times its rest mass. The 
rest mass of the nucleus (the proton) could be considered as the lowest limit of 
a macroscopic body, being in this case the electron the microscopic particle. 

If the principle of this phenomenon, for the formation of UHECR in colli- 
mated jets, is realistic, we have to expect the total mass of the ingoing body to 
be "macroscopic" as compared to (i.e. infinitely larger than) the total mass of 
the ejected particles, and to be of the same order of magnitude than the mass 
falling in the BH after each Penrose type decay inside the ergosphere. 

Then, in terms of "rest mass" (i.e. of the amount of matter), the efficiency of 
the Penrose process is very (infinitely) weak, meaning that a very weak part of 
the ingoing mass is ejected to the outside. But in terms of energy (a conserved 
quantity), the Penrose process is very (infinitely) efficient, meaning that the 
energy of the ejected particles is much (infinitely) larger than the energy of 
these same incident microscopic particles (when they were linked yet to the 
incident macroscopic body). 

We can make an analogy with a " superelastic" shock, as sometimes occurs in 
nuclear physics, in which the total kinetic energy of the final outgoing particles 
is greater than the kinetic energy of the initial ingoing particles, the difference 
being acquired to the detriment of a part of the ingoing (rest) mass, called 
"mass defect" (as opposed to an "elastic shock" , for which the kinetic energy is 
conserved, or to an "inelastic shock" for which the final kinetic energy is lower 
than the initial one) . In our case, the lost mass falling into the BH is equivalent 
to the " mass defect" , and is found again in the form of kinetic energy of the 
ejected mass. Besides, the efficiency of a collimated jet formation (in terms 
of mass as well as of energy) is apparently even more weakened by the fact 
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that, on the set of ejected particles (in all the directions from the ergosphere, 
by this process), very few are ejected in a direction which is asymptotically 
parallel to the z-axis [23] ■ However, the only possible collimation being in this 
direction, the coherent set of these particles gives them a privileged observable 
character, which does not exist for the other particles isolatedly scattered in all 
the directions. 

Recent results of the Pierre Auger Observatory [33] seem to show a correla- 
tion between the UHECR (above 57EeV) and the nearby (< 71 Mpc)) AGNs. 
It would be interesting to examine more precisely the correlations with the di- 
rections (radioloud, quasars or blazars) of the jets and with the proximity of 
their sources. 

All the results, up to the section 4 (in particular the numerical values, which 
concern adimensional quantities), are strict consequences of the structure of the 
Kerr metric where we solely fixed the BH spin, a — M/2, and we made the 
assumption of a double real root Y. 

A first expected position of the asymptote corresponding to an "infinite" 
energy in our model is p\ ~ 10p e = 5M = 2.5 * rg. The second asymptote pre- 
dicted by our model for a very high energy is p\ = (0.693199/2)M ~ 0.3466M, 
i.e. at about 1/6 of the Schwarzschild radius r$ — 2M. This last case is the 
only compatible with the limitations imposed by the characteristics of the sys- 
tem of geodesies equations, as seen in section 6. That means that the presently 
described jet is a very thin (narrow) jet. 

Besides, there will be a strong concentration of the most energetic part of 
the beam in the close vicinity of these two asymptotes, with, on one side, an 
abrupt (steep) decreasing and, on the other side, a smoother decreasing, which 
indicate a radial energetic structure of the jet, as observed [13] . 

The rarity of the currently detected UHECR [3] could be explained in our 
model by the rarity of UHE particles in the very narrow beams near the two 
asymptotes shaping the jet, i.e. directly observable in the very precise ax- 
ial direction only. Even though the energies have to be limited (< 50EeV) 
by the Greisen-Zatsepin-Kuzmin effect, firstly primary particles with energies 
< 50EeV could come from nearby (< 70 MPc) AGN to the earth (principally 
neutrinos), and secondly, for more energetic and more distant primary particles, 
a particular mechanism, for instance with secondary particles, can be consid- 
ered p~2] • Anyway, the recent observations of multi-TeV photons from distant 
blazars which do not display the expected spectrum (with suppressions due to 
the interaction with the cosmic microwave background) require an explanation 

®- 

Let us briefly recall that we found these results under some restrictive condi- 
tions, taken into account in section 3, especially the assumption of the existence 
of a double root Y, and the choice of the "middle" value a — M/2. We can thus 
hope that by "relaxing" these assumptions, other more general results could 
emerge. 

For instance, relaxing the assumption of a double root Y, while keeping 
the same value M/2 of the parameter a, could open the possibility of other 
solutions with pi > p e , i.e. thicker jets, for admissible high energy jets. Work is 
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in progress in that way and our first results are encouraging. Another example 
could be to study, in this model, the role of the BH spin a by trying (say) 
higher values of this parameter. Let us note, by the way, that there are few 
observational results concerning the possible values of a [37J [501 13 IS] > an d its 
possible correlation with the length and power of the jet, while there are more 
numerous observations concerning the values of M. 

Our results can also easily be extended to particles other than electrons, for 
example to protons or neutrinos. This does not change the "geometry" that we 
obtained, i.e. the positions p\ of the jets, but their energy only. 

For a proton {-/Si ~ lGeV) the maximal energy we can here numerically 
calculate (but which is theoretically as large as we want) is about E ~ 5.6 * 
\Q 25 eV — 5.6*10 7 EeV, which largely includes the highest energies of the current 
observed UHECR [3 Ql [19]. 

Finally, let us also recall that our model does not require magnetic fields, 
which allows us to discard some problems related to their strength [FT] , necessary 
to obtain such huge energies, and permits us to consider neutral particles as well. 
For example, a neutrino, which mass is assumed to be y/Si — 0.33ey [35], for 
the precedent evaluation, would reach the energy E ~ 2*10~ 2 EeV, which seems 
to be an acceptable value [HI HH [8J . 

In our model, the only role played by a magnetic field is the role of the 
induced magnetic field which tends to stabilize the jet (e.g. [2] [21] ). The 
inner part of the parallel jet, if composed for instance of electrons only, creates 
a magnetic field, which tends to stabilize the collimation for its outer part 
(which is the most energetic part in the precedent example p\ — 0.3466M). 
The strength of the magnetic field depends on the relative part of the charged 
particles in the jet. The radial structure of the jet, in our model, supported 
by direct observations [15], would require a stability study. However, a first 
rough evaluation shows that the ratio \[S\jEc — yfix/E remains always <C 1 
inside the UR jet, which let us think that any Kelvin-Hclmholtz instability is 
negligible, in accordance with the results of extensive studies in 2D [25] or in 
3D [1FJ. 

Recent articles discussed the possibility to generate high energy particles by 
collisions near a BH [UJ and evoked the interest to consider the Penrose process 
[7]. Our approach can be seen as a contribution to this debate. 
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